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ABSTRACT 



The integrated bolometric effective surface brightness distributions of starbursts 
are investigated for samples observed in 1. the rest frame ultraviolet (UV), 2. the far- 
infrared and Ha, and 3. 21cm radio continuum emission. For the UV sample we exploit 
a tight empirical relationship between UV reddening and extinction to recover the bolo- 
metric flux. Parameterizing the upper limit by the 90th percentile of the distribution, 
we find a mean S'e^go = 2.0 x 10^^ Lq kpc~^ for the three samples, with a factor of three 
difference between the samples. This is consistent with what is expected from the 
calibration uncertainties alone. We find little variation in 6'e,9o with effective radii for 
i?e ~ 0.1 — 10 kpc, and little evolution out to redshifts 2; ~ 3. The lack of a strong 
dependence of 5'e,9o on wavelength, and its consistency with the pressure measured in 
strong galactic winds, argue that it corresponds to a global star formation intensity limit 
(Se,9o ~ 45 kpc~^ yr"-*^) rather than being an opacity effect. There are several im- 
portant implications of these results: 1. There is a robust physical mechanism limiting 
starburst intensity. We note that starbursts have consistent with the expectations of 
gravitational instability models applied to the solid body rotation portion of galaxies. 
2. Elliptical galaxies and spiral bulges can plausibly be built with maximum intensity 
bursts, while normal spiral disks can not. 3. The UV extinction of high- 2; galaxies is 
significant, implying that star formation in the early universe is moderately obscured. 
After correcting for extinction, the observed metal production rate at 2; ~ 3 agrees well 
with independent estimates made for the epoch of elliptical galaxy formation. 

Subject headings: galaxies: starburst - ultraviolet: galaxies - infrared: galaxies - radio 
continuum: galaxies - early universe 
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1. Introduction 

Starbursts are regions of intense massive star 
formation that can totally dominate a galaxy's 
integrated spectrum. They range in size from gi- 
ant extragalactic Hll regions of scale size ~ 0.1 
kpc (e.g. NGC 604 in M 33) to global bursts (of- 
ten in merging systems) many kpc across that 
can cover the entire face of the system (e.g. 
NGC 4038/4039, NGC 4449). There has been 
much interest in the starburst phenomenon over 
the past decade or so, yet there is much more to 
learn. For example we do not know how star- 
bursts are regulated: what turns them on and 
off, and what keeps them going. Also unknown 
is how the properties of starbursts vary with size 
or luminosity. We address some of these issues 
in this paper. 

While starbursts are worthy of investigation 
in their own right, they are even more important 
when placed in the broader context of contempo- 
rary extragalactic astrophysics. The cosmologi- 
cal relevance of starbursts has been underscored 
by the recent discovery of: the existence of a pop- 
ulation of high-redshift (z > 2) UV-bright field 
galaxies (cf. Steidel et al. 1996a,b; Lowenthal 
et al. 1997). The sheer number density of these 
galaxies implies that they almost certainly rep- 
resent precursors of typical present-day galaxies 
in an early actively-star-forming phase. This dis- 
covery moves the study of the star-forming his- 
tory of the universe into the arena of direct obser- 
vations (cf. Madau et al. 1996), and gives added 
impetus to the quest to understand local star- 
bursts. In particular, a thorough understanding 
of how to exploit the diagnostic power of the rest- 
frame ultraviolet (UV) properties of local star- 
bursts will give astronomers powerful tools with 
which to study star-formation in the early uni- 
verse. 

Accordingly, we have been using the HST to 
determine the basic UV structure and morphol- 
ogy of starbursts. In our first paper (Meurer 
et al. 1995; hereafter M95) we analyzed in de- 



tail UV HST images of nine starburst systems 
with distances^ D < 75 Mpc. Here, we turn our 
gaze outward and consider the UV structure of 
starbursts from z = to z > 3, in an attempt to 
document and understand the sytematic proper- 
ties of starbursts. We also broaden our net by 
considering the structural properties of infrared- 
selected "dusty" starbursts. 

M95 showed that nearby starbursts are irreg- 
ular structures consisting of diffusely distributed 
light interspersed with prominent compact (radii 
< 10 pc) star clusters. Numerous other HST 
studies have also commented on the presence of 
luminous young star clusters (including Holtz- 
man et al., 1992; Whitmore et al, 1993; Conti & 
Vacca, 1994; Hunter et al, 1994; Bower & Wil- 
son, 1995; O'Connell et al, 1995; Whitmore k 
Schweizer, 1995; Holtzman et al., 1996; Maoz, 
et al. 1996a, b; Schweizer et al, 1996; Watson 
et al, 1996). While the clusters are certainly 
the most striking aspect of the images, it is the 
diffuse light that dominates. M95 find on average 
80% of the UV light is diffusely distributed, and 
20% comes from clusters. Maoz et al. (1996a) 
find similar fractional UV contributions of clus- 
ters in five starburst ring galaxies. M95 ar- 
gue that there are two modes of star formation 
in starbursts: prominent cluster formation, and 
dominant diffusely distributed star formation. 

M95 also noted that the effective surface bright- 
nesses of most of the starbursts in their (small) 
sample span a narrow range of values (disper- 
sion of 0.4 mag arcsec"^ for 8/11 starburst re- 
gions). Since starbursts produce the most in- 
tense UV emission observed in galaxies, this re- 
sult implies that there is an upper limit to the 
UV surface brightness of galaxies. Lehnert &; 
Heckman (1996) found a similar limit to the far- 
infrared surface brightness of far-infrared galax- 
ies (FIRGs). In both samples the illuminating 
source can be traced to high-mass stars. The 

^We convert all results to Ho — 50kms^^ Mpc^^,Qo ~ 0, 
except where noted. 
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UV light studied by M95 is dominated by the 
hght from stars with mass ~ 20 Mq. The far- 
infrared emission studied by Lehnert & Heckman 
is thought to result from dust absorbing UV - op- 
tical photons (which are predominantly produced 
by high-mass stars) being heated and reemitting 
the radiation thermally at infrared wavelengths. 
Thus FIRGs are thought to be obscured star- 
bursts. These surface brightness limits suggest 
a limit to the intensity of high-mass star forma- 
tion, which implies that there is a mechanism 
that limits the global star formation intensity of 
galaxies. 

Here we reexamine the issue of starburst inten- 
sities. Since selection effects are strongest at low 
surface brightness, we are primarily interested in 
the intensity maximum. Our aims are to deter- 
mine what the starburst intensity limit is, inves- 
tigate whether it varies with certain parameters, 
and consider what causes the limit. By com- 
paring starbursts selected and observed at differ- 
ent wavelengths we address the issue of opacity 
and the difference between dusty and relatively 
dust free starbursts. By comparing starbursts 
at different redshifts we are able to consider the 
evolution of starburst intensities, and the resul- 
tant cosmological implications. Our method is 
to compile observations of a diverse set of star- 
burst samples from the literature, combined with 
(a dash of) new data. After defining some rele- 
vant quantities in §^ we compile three samples 
of starburst data: I. a sample of starbursts ob- 
served in the rest-frame vacuum ultraviolet (§^; 
II. a sample observed in the far-infrared and Ha 
(§^; and III. a sample observed at radio wave- 
lengths (§|5|). In §|6|we discuss the selection effects 
and measurement biases, and estimate the level 
of agreement we expect between intrinsically sim- 
ilar samples. In we discuss the interpretation 
of the results and their implications, and §|8| sum- 
marizes the conclusions. 



2. Definition of Quantities 

To quantify the intensity of high-mass star for- 
mation we measure the surface brightness at a 
wavelength that traces high-mass stars, within 
an aperture that encompasses half of the total 
emission from the starburst. Thus we employ the 
effective, or half light, radius Re, and effective 
surface brightness enclosed within R^. These 
quantities are best measured with a curve of 
growth extracted from a background-subtracted 
image of the starburst. Relative to the total 
luminosity L, the effective surface brightness is 
given by: 



Note that the definition of surface brightness 
used by Lehnert & Heckman (1996) is a factor of 
two higher than our definition. Their definition 
is equivalent to saying that all of the starburst is 
located within the half light radius of the tracer. 
Here we assume a one-to-one spatial relationship 
between the high-mass stars and their trace emis- 
sion. One complication is that the apertures used 
to measure Re are often elliptical, not circular, in 
shape. Our approach is to take Re as the geomet- 
ric equivalent radius of the aperture: 

Re = ae^{h/a)e (2) 

where is the semimajor axis length of the el- 
liptical aperture enclosing half the the total flux, 
and (b/a)e is the corresponding axial ratio at Og. 

We emphasize that Re and Se are global quan- 
tities; that is they are associated with the total 
light of the starburst. This includes both clus- 
ters and diffuse light. Furthermore the area over 
which the Se is measured scales with the distribu- 
tion of the light. For moderate to high redshift 
galaxies full cosmological corrections have been 
applied to the Re and L measurements assum- 
ing an open universe with Hq = 50 km s~^ Mpc~^ 
and Qo = 0. 

The salient properties of starbursts are es- 
sentially defined by ionizing stars. These have 
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masses ^ 20 Mq, and hence lifetimes <^ 10 Myr. 
Ideally we would like to determine how these 
stars are distributed. However, the size of the 
ionizing source distribution is not directly ob- 
servable since either very little ionizing emission 
escapes from starbursts, or it is absorbed by the 
local ISM/IGM before reaching the earth (Lei- 
therer et al, 1996). Therefore Re of some other 
tracer of the ionizing population is employed: 
the vacuum- UV (sample I), Ha (sample II), and 
21cm continuum radio emission (sample HI). To 
compare results we convert all flux related mea- 
surements to bolometric quantities. The algo- 
rithms we use are given below. 

To characterize the Se distribution of a sam- 
ple we consider the sample median and 90th 
percentiles, which we denote as 5*6,50 and 5e,9o, 
respectively. The samples were compiled from 
a hodge-podge of studies with selection criteria 
and measurement techniques varying from study 
to study. In addition, there is a well known 
bias when observing extended objects to observe 
those with the highest surface brightness (e.g. 
Disney 1976; McGaugh et al. 1995), so the lower 
percentiles of the Se distributions are likely to 
be rather incomplete. This is especially true for 
starbursts which are often recognized by their 
high surface brightness. Therefore, /Se.so and 
Sg^QO should not be interpreted too literally. They 
actually represent upper limits to the true me- 
dian and 90th percentiles of the Se distributions 
of star forming galaxies. Here 5e,5o should be 
interpreted as the typical surface brightness se- 
lected in a sample, and Se,9o should be inter- 
preted as close to the highest surface brightness 
found by a given technique. In §^ we discuss the 
selection effects and biases of the samples and 
estimate the level of agreement we expect. 

3. Ultraviolet sample 
3.1. Ultraviolet extinction 

There are many appealing reasons to observe 
starbursts in the vacuum- UV (see M95). How- 



ever, there is one serious hindrance to interpret- 
ing the results: dust. It efficiently extincts and 
scatters UV radiation. When dust absorbs UV 
and optical photons it heats up and reemits the 
radiation thermally in the far-infrared. The re- 
distribution of radiant energy must be well mod- 
eled and corrected for in order to determine the 
true bolometric luminosity of a starburst. 

M95 showed that it is possible to do just that 
using only the UV properties of starbursts. This 
is illustrated in Fig. || which is adapted from 
Fig. 6 of M95. It shows the ratio of far- infrared 
(FIR) flux FpiR to UV flux F220 compared to the 
spectral slope /3 which is a measure of the ul- 
traviolet color, for a sample of UV-selected star- 
bursts. The definitions of these quantities are 
given below (§4.1, § |3.2| ) and in M95. In this plot 



the y axis quantifies the redistribution of spec- 
tral energy from the UV to the FIR. A strong 
relationship is apparent in this diagram - in the 
sense that, as relatively more light is emitted in 
the FIR, the starbursts become redder. Note 
that the hatched region in this plot shows the ex- 
pected range of colors for naked ionizing popula- 
tions. All the galaxies in this plot display strong 
recombination emission spectra, hence their in- 
trinsic colors should be that of an ionizing popu- 
lation. Figure || then indicates that reddening is 
positively correlated with dust extinction. Such 
a correlation is a major prediction of simple fore- 
ground screen models for dust extinction. The 
dashed line in Fig. |l| shows the expected rela- 
tionship for a starburst having /? = —2.5 and 
foreground screen dust having the Calzetti et al. 
(1994; hereafter C94) "extinction law"0. Details 
of the model can be found in M95. 

Figure |l| shows that a foreground screen dust 
geometry models well the redistribution of spec- 



''We use the term "extinction law" rather loosely. The 
C94 law is better referred to as an attenuation law since 
it is based on observations which recover most or all of 
the UV flux, hence extinction is partially compensated for 
by foreward scattering, thus producing a relatively grey 
attenuation law. 
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tral energy from the UV to the FIR. Similar re- 
sults have been noted by Calzetti et al. (1996) 
and Lehnert & Heckman (1996). However, dust 
in the vicinity of starbursts does not necessarily 
have such a simple distribution in all cases. The 
rms in -Ffir/^220 about the model line of ~ 0.4 
dex is probably larger than the errors, and may 
be indicative of the non-uniform nature of the 
dust distribution (e.g. dust lanes). In addition, 
Fig. H plots only UV-selected starbursts, which 
preferentially will have little extinction. 

The point is that the tight empirical relation 
between UV color and infrared to UV flux ratio 
shown in Fig. || gives us a method to remove the 
effects of dust extinction for UV-selected star- 
bursts. We adopt the C94 curve in Fig. |l| be- 
cause it is physically plausible and goes through 
the data fairly well. Whether or not the details 
of this model are absolutely correct is irrelevant, 
as long as this empirical relationship holds for 
all UV-selected starbursts. Other plausible fits 
through the data could also be used, and they 
will yield the same UV correction factor as we 
derive to within 0.4 dex. 

3.2. Calculation of intrinsic quantities 

In all cases measurements of L, R^, and thus 
Se have been made from HST images obtained in 
the rest frame UV. In order to make extinction 
and k corrections it is crucial to measure the UV 
color, or equivalently the UV spectral slope, since 
the observed UV continuum spectra of starbursts 
are well fit by a power law of index /3: 

fx oc (3) 

(C94, M95). Here f\ is the spectral flux density 
per wavelength interval. In order to compare the 
subsamples we reference all observations to the 
observations of M95. These employed the Faint 
Object Camera with F220W filter, which has cen- 
tral wavelength Ac = 2320 A. The total flux is 
given as F220 = Ac/220 • The methods used to es- 
timate f3 vary from subsample to subsample, and 
are described below. 



The first step in converting the UV to a bolo- 
metric flux is to apply (where appropriate) the k 
correction: 

,(220) = ^ = |I1±4^|''. (4) 

JXc/{l + z) [ Ac J 

Here Ac corresponds to the central wavelength 
of the filters used below. After correcting for 
Galactic foreground extinction, the UV fluxes 
were then corrected for intrinsic extinction using 
the C94 model line in Fig. |^ as discussed above. 
This model yields the intrinsic extinction ^int fol- 
lowing the algorithm in M95. For this model a 
change in spectral slope A/3 = 0.1 corresponds 
to A^int ~ 0.2 mag for small (3 — Pq. 

We use the population models of Leitherer & 
Heckman (1995; hereafter LH95) to convert from 
intrinsic F220W flux/luminosity to bolometric 
flux. We assume that starbursts are well rep- 
resented by a solar metallicity 10 Myr duration 
constant star formation rate population having a 
Salpeter (1955) initial mass function (IMF) span- 
ning the mass range of 0.1 to 100 Mq. The 
10 Myr duration is meant to represent the typi- 
cal crossing time of a starburst as estimated by 
M95. Although longer durations are likely in 
the largest starbursts, the UV properties will be 
dominated by the stars less than 10 Myr old. In 
any case the choice of a particular age has only 
a minor effect on the determination of bolomet- 
ric quantities. This is demonstrated in Fig. ^ 
which shows the evolution of i22o/-^boi at so- 
lar metallicity for both instantaneous burst and 
constant star formation rate populations for this 
IMF slope (see M95 and LH95 to see the evo- 
lution of the numerator and denominator sepa- 
rately). For our adopted population model, 



Figure |2| shows that for this IMF slope the adopted 
bolometric correction is accurate to about 20% or 
better for instantaneous burst populations less 
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than 10 Myr old (i.e., ionizing bursts) and con- 
stant star formation rate populations over the 
full age range shown. For our adopted starburst 
population parameters the mass to light ratio is 



M 

Lhoi 



2.2 X 10 



-3 



(6) 



and the relationship between star formation rate 
M. and Lboi is 



L 



bol 



M 



4.5 X 10^0 ITWoyr 



-1 ■ 



(7) 



where the luminosity is referenced to Lq = 
3.83 X lO^^ergs"^, the bolometric luminosity of 
the sun (Allen, 1973). The conversions of light 
to mass (eq. |6|) and star formation rate (eq. |^ 
are much more dependent on the IMF parame- 
ters and adopted star formation history than the 
bolometric correction (eq. |5|), hence our mass and 
star formation rate estimates should be consid- 
ered indicative only. 

3.3. UV Subsamples 

(La) Local UV-selected starbursts are repre- 
sented by the M95 dataset, where the results are 
derived from HST Faint Object Camera (FOC) 
images obtained with the F220W filter. Total 
fluxes and effective radii were determined for the 
eleven starburst regions observed m nine galax- 
ies (there are three detached regions in their 
NGC 3690 image). The sample is at sufficiently 
small redshift that k corrections are unneces- 
sary. The Galactic foreground extinction is taken 
from Burstein & Heiles (1984), and removed from 
the UV following the extinction curve of Seaton 
(1979) using the algorithm given by M95. The f3 
values were measured from International Ultra- 
violet Explorer spectra, and taken directly from 
M95. 

(Lb) We have obtained new data on two star- 
burst galaxies at z ^ 0.4, QNY1:32 and SGP1:10, 
discovered by Boyle et al. (1990). The two galax- 
ies were observed with the HST using the FOC 
and the F342W filter (Ac = 3403A), which has a 



central wavelength in the rest frame of the galax- 
ies {Xc/{l + z) ^ 2380 A), very similar to that of 
the F220W filter. The k corrections (eq. g), al- 
though small (< 0.1 mag), were applied. Both 
galaxies consist of three closely spaced knots 
aligned nearly linearly. The 512 x 512 images 
were corrected for non-linearity of the FOC and 
the fluxes and effective radii were obtained from 
curves of growth using concentric elliptical aper- 
tures. The apertures were centered on the cen- 
tral knot and had orientations and axial ratios 
set to match the outer isophotes of the galaxy. 
The basic properties of the galaxies are listed in 
Table |. 

Galactic extinction, listed as ^cai in table ||, 
was estimated from the foreground Hi column 
density (Stark et al. 1992) using the Mathis 
(1990) extinction curve. These galaxies do not 
have (3 or UV color measurements, so we have to 
estimate ^int by some other manner. We have 
obtained flux calibrated optical spectra for both 
sources, which confirms that these are indeed 
starbursts (strong narrow emission lines). The 
SGP1:10 spectrum was obtained by us with the 
KPNO 4m telescope. Alexei Filippenko kindly 
obtained the QNY1:32 spectrum for us using 
the Lick 3m telescope. The ratio Fn-y/Fnp = 
0.46 lb 0.11 for SGP1:10, which is consistent with 
case B recombination and no intrinsic extinction, 
hence we assume ^int = 0. Only one Balmer 
line (H/3) was well detected in the spectrum of 
QNY1:32. We assume it has the median redden- 
ing (/? — Pq) = 1.4 of the rest of the UV sample, 
hence ^int = 2.2 mag. 

(I.c-e) Very high redshift galaxies are repre- 
sented by three samples. One was selected from 
ground based observations, and the other two 
were from the Hubble Deep Field (HDF) ob- 
servations (Williams et al, 1996). In all three 
cases the galaxies were selected as probable Ly- 
man break systems from their broad band col- 
ors using the technique pioneered by C. Stei- 
del and collaborators. We consider only the 
galaxies that have been spectroscopically con- 
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firmed, which has been done with the Keck tele- 
scope in all cases. The confirmation spectra fre- 
quently showC IV and Si IV absorption features 
which are commonly observed in starburst galax- 
ies. These arise in the winds and photospheres 
of high-mass stars and/or a highly ionized in- 
terstellar medium. The crucial size measure- 
ments come from HST imagery with WFPC2. 
The spectral slope /?, is measured directly from 
published broad band colors in the AB system 
(mAB = —48.6 — 2.5 log(/;y), where /j^ is the spec- 
tral flux density per frequency interval in units 
of ergcm"^ s^^ Hz^^). The redshifts {z ~ 3) of 
the sources put the V to I band observations 
squarely in the rest frame UV. Details on the in- 
dividual high-z subsamples are as follows: 

(I.c) The original selection for this subsam- 
ple is from the ground based photometry of Stei- 
del k Hamilton (1992) and Steidel et al. (1995). 
Spectroscopic redshifts are given by Steidel et al. 
(1996a), and HST imagery is reported by Gi- 
avalisco et al. (1996). There are seven spectro- 
scopically confirmed objects in the sample. Gi- 
avalisco et al. show that two of the seven galax- 
ies are double. Here we consider the components 
separately assuming they have a common red- 
shift, so the total size of this subsample is nine. 
We derive j3 from the {Q — 7^)ab colors given by 
Steidel et al. (1996a): 

/3 = 2.55(e? - 11) AB - 2, (8) 

We take = -^1/2 from Table 2 of Giavalisco 
et al. (1996). 

(I.d) The first HDF sample is that of Steidel 
et al. (1996b). Five confirming redshifts were 
obtained, one of which corresponds to a double 
source; thus the total subsample size is six. Ef- 
fective radii Re = R1/2 were taken from their Ta- 
ble 2. Rest-frame UV fluxes were derived from 
the 7^ band (combined F606W and F814W hght) 
magnitudes in their Table 1. Finally, P was es- 
timated from the (rn-Feoew — '71f814w)ab colors 
given in the HDF version 2 catalog (Williams 



et al, 1996): 

P = 3.23(mF606W - m'F814w)AB - 2. (9) 

We used (mpeoew — "T'F814w)ab in preference to 
(mF45ow — '^)ab listed by Steidel et al., because 
the blue band in the latter may be affected by 
opacity from the Lya forest (Madau et al, 1996). 

(I.e) The other HDF subsample has confir- 
mation spectra from the DEEP ("Deep Extra- 
galactic Evolutionary Probe") program (Lowen- 
thal et al. 1997). A total of 11 high-z detections 
were made for sources not in sample I.d. Four of 
these belong to double sources. Hence the total 
subsample size is 15. Effective radii were taken 
from their Table 2 (listed as and /? was 

derived from the (mFeoew — ™'F814w)ab colors 
given in their Table 1 using eq. |9[ 

3.4. UV results 

Figure ^ shows Lboi and Se clS 3j function Rp 
of the UV-selected samples. The 5*6,50 and Se^go 
surface brightness levels of the combined sam- 
ple are plotted as dashed and dotted lines re- 
spectively. Figure ^ shows the surface bright- 
ness distribution of the combined UV sample 
in the top panel. The combined sample has 
Se,9o = 2.0 X lO"L0kpc-2. Table | tabulates 
log(5e,5o) and log(5e,9o) levels of the combined 
UV sample, and its various subsamples. 

It is immediately apparent from Fig. ^ that Se 
shows little or no dependence on Re over about 
two orders of magnitude in size; hence there is 
no dependence on Lboi over about four orders of 
magnitude in luminosity. Table |l] indicates that 
/Se.go is the same within 0.1 dex for the entire 
sample and the combined subsamples I.c-e con- 
taining just high-z galaxies. This is quite small 
compared to the cosmological dimming (which 
has been removed) of log(z -|- 1)'* = 2.4 dex at 
z = 3. 

All of the starbursts in this sample are well 
resolved; there are no upper limits on size, even 
though the samples were not biased against point 



8 



sources. The range of angular sizes of the galax- 
ies in the samples spans nearly two orders of mag- 
nitude similar to the range of absolute sizes. We 
checked how our results are affected by whether 
detached substructures are counted separately or 
not by considering the two galaxies in subsample 
Lb. When their knots are considered separately 
the brightest in each has about 0.25 dex higher 
than the galaxy considered as a single unit. How- 
ever, the average Se of the individual knots is 
close to the Se value of the galaxy as a whole. 
This indicates that the detached substructures 
have a high covering factor. We conclude that 
the Se measurements are robust to within ~ 0.3 
dex with respect to resolution effects. 

The same can not be said about the embedded 
structure. In particular the star clusters within 
starbursts have much higher Se values than star- 
bursts clS Si whole. Figure ^ shows Lboi and Se 
of the star clusters with size estimates from M95 
(their table 10). The bolometric corrections were 
estimated using eq. ^. The ^e.so and Se,9o lines 
from Fig. ^ are plotted for comparison. All of 
these clusters have surface brightnesses above the 
<S'e,5o and most are above the 5*6,90 level derived 
from integrated starbursts. This is partially a 
selection effect since the clusters are all embed- 
ded within the dominant diffuse background of 
a starburst. In this figure there are many up- 
per limits to Re, and thus lower limits to Se- In 
addition the Re estimates for the clusters in the 
more distant galaxies are probably inflated (and 
Se underestimated) as noted by M95. The most 
intense emission from a resolved cluster comes 
from NGC1705-1 with = 4.7 x 10^^ Lq kpc'^. 
This is over two orders of magnitude more intense 
than 5e,9o of the UV sample, yet NGCl 705-1 has 
probably already faded by a factor of ~ 6 from 
when it was an ionizing cluster (M95, Meurer 
et ai, 1992). Figure |5| demonstrates that high Se 
sources can be detected in the UV, they just do 
not have large scale sizes, nor do they dominate 
the integrated light of starbursts. 



4. Far-infrared - Ha sample 
4.1. Subsamples, method 

The two subsamples of FIRGs are (II. f) the 
sample of Armus, Heckman &: Miley (1990); and 
(II. g) the sample of Lehnert & Heckman (1995, 
1996). They were drawn from the IRAS catalog 
using the same FIR flux limits used to create the 
IRAS Bright Galaxy Survey (BGS; Soifer et ai, 
1987, 1989). In addition, these samples have 
an infrared color selection, /eo/Zioo ^ 0.77,0.4, 
respectively, where /so, /loo are the 60/im and 
lOO/um flux densities from the IRAS Point Source 
Catalog (1988). In order to minimize contamina- 
tion of the samples by galaxies dominated by ac- 
tive galactic nuclei (AGN) we excluded galaxies 
with a LINER or Seyfert spectral classification 
in Veilleux et al. (1995). 

The FIR luminosity is derived from the far- 
infrared flux given by 

Ffir = 1.26 X 10-1^2.58/60 + /100) erg cm-2 s'S 

(10) 

and the units for /go, /loo are Jy (Helou et al., 
1985). Almost all of the galaxies in these samples 
have have Ffir > -^optical- Hence to a good ap- 
proximation the bolometric luminosity is equiv- 
alent to that produced in the FIR. We take 

^boi = c(/6o//ioo)-Pfir, (11) 

where c(/6o//ioo) is the bolometric correction to 
-^FiR to obtain the total FIR luminosity (wave- 
length range 1-1000 ^m). The correction as- 
sumes a single dust temperature and emissivity 
(X and is tabulated by Helou et al. (1988). 
For subsample II. f we assume a single average 
c(/6o//ioo) = 1-5, on the basis of the /eo/Aoo 
color range of the sample. For subsample II. g we 
use the galaxy by galaxy c(/6o//ioo) corrected 
fluxes given by Lehnert and Heckman (1996). 
Since the Lboi measurements are made in the 
FIR, no extinction corrections are necessary. For 
the size of the starburst we take Re = i?e(Ha), 
the Ha effective radius. 
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The Ha images used to estimate Re were ob- 
tained from the ground in FWHM seeing typ- 
ically 1-2". Therefore we considered measure- 
ments with Re < 1.0" to be upper limits. There 
is only one galaxy in the two subsamples with 
such an upper limit: 01217+0122 with Re < 
0.86" from subsample II. f. Its actual Sg value 
does not significantly affect our results because 
its limiting Se > 10^^ Lq kpc"^ is already greater 
than 5e,9o of the combined FIR/Ha sample (see 
below) . 

4.2. FIR/Ha results 

Figure ^ shows Lboi and Se function of 
Re for sample II. In comparison to the UV sam- 
ple, shown in Fig. |3|, the FIR/Ha sample looks 
less like a constant surface brightness correla- 
tion. The FIR/Hq sample is richer in low sur- 
face brightness galaxies than the UV sample, as 
can be seen from the surface brightness distribu- 
tion shown in the second panel of Fig. ^ How- 
ever, there is little apparent variation in the up- 
per envelope of points in this sample. This 
illustrates nicely that we are dealing with a sur- 
face brightness limit and not a constant surface 
brightness correlation. 

The combined FIR/Ha sample has ^e^go = 
0.83 X 10"LQkpc-2. Table [l| tabulates the 
log(S'e,5o) and log(S'e,9o) values for the two sub- 
samples and the combined FIR/Hq sample. These 
values for the combined sample are shown as 
dashed and dotted lines respectively in Fig. ^ 

5. Radio sample 
5.1. Sample, method 

In the radio we use the sample of Condon et al. 
(1990; hereafter C90) which was originally drawn 
from the IRAS BGS. The C90 sample includes 
over 300 IRAS galaxies of all types observed at 
21cm with the Very Large Array. The editing 
steps we apply, which are discussed below, reduce 
this to a final sample of 38 sources. 



To determine the bolometric flux we use the 
FIR to radio correlation. This well-known cor- 
relation (see Helou, 1991, for a review) between 
fluxes in the two wavelength regimes holds for 
over four decades in luminosity with an intrin- 
sic scatter of only 0.2 dex making it the tightest 
correlation in global fluxes of galaxies. Although 
the physical explanation of the correlation is only 
now being addressed (e.g. Lisenfeld et ai, 1996), 
it seems to work remarkably well for star forma- 
tion in all sorts of environments including normal 
disk galaxies, starburst systems, and even resid- 
ual star formation in elliptical galaxies. The form 
of the FIR - radio correlation is: 



-^FIR 



3.75 X 1012 Hz 



10V./(21cm) 



(12) 



where q ~ 2.35 for starburst and normal galaxies 
(Sanders & Mirabel, 1996). We adopt this rela- 
tionship to estimate -Ffir and then obtain Fboi 



by bolometric correcting FpiR using eq. |11| and 
adopting c(/6o//ioo) = 1-4. This is appropriate 
for the /eo/Zioo ratios of our final radio sample 
and corresponds to a dust temperature of T = 45 
K to 70 K, and a i^"^ emissivity (Helou et al, 
1988). 

COO report 2D-Gaussian fits to the radio im- 
ages, often for a variety of spatial resolutions and 
using multiple components. We derive Re from 
the angular effective radius taken to be 



A/ Cm 



(13) 



where 9m and 9m are the major and minor axis 
FWHM of the 2D fits. For a true circularly sym- 
metric Gaussian distribution 9e = 9m/2 = 9m f^- 
Condon et al. (1991) present higher resolution 
A = 3.6 cm sizes of ultr aluminous FIRGs. We 
chose not to use this study because for most of 
the galaxies in common to our final sample, the 
components they resolve only comprise a small 
fraction of the 3.6 cm flux, or comparisons of the 
21 cm and 3.6 cm maps suggest that much dif- 
fuse emission is missing from the 3.6 cm maps. 
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The few remaining galaxies have 9e that agree to 
within a factor of ~ 2 at 3.6 and 21 cm. 

Four editing passes were done to the sample 
in order to use well resolved observations and 
isolate the primary contributor to systems most 
likely to be starbursts (i.e. not AGN). Firstly, for 
a given spatial resolution only components with 
/,y(21cm)//,y (21cm, total) > 0.5 were selected so 
as to reject minor contributors to the total flux. 
Secondly, for systems observed with with more 
than one VLA configuration (i.e. at different res- 
olutions) we selected the best resolved observa- 
tions by taking those with the largest Qu^ml 
[Beam area]. Usually this meant selecting the 
highest resolution observations. Thirdly, we re- 
ject galaxies that do not obey the radio-FIR cor- 
relation; from the IRAS and 21 cm fluxes tabu- 



lated by 090, we calculate q in eq. |I2| and reject 
sources with q < 1.9 as likely radio loud AGN 
(Sanders & Mirabel, 1996). Finally, in order to 
exclude other known AGN, we further selected 
only the systems with Hll type spectra in the nu- 
clear regions according to Veilleux et al. (1995). 

5.2. Radio results 

The resulting sample has a relatively narrow 
1.5 dex luminosity range centered on Lboi = 
2 X 10^^ Lq. This is primarily due to the fi- 
nal editing step. The BGS subsample used by 
Veilleux et al. (1995) is weighted towards galax- 
ies with Lboi ^ 10^^ (Kim et al. 1995), hence 
low luminosity galaxies were discarded. Rejec- 
tion of AGN preferentially removes the high- 
est luminosity IRAS galaxies (Veilleux et al, 
1995; Sanders & Mirabel, 1996). The surface 
brightness distribution of this sample is shown 
in the third panel of Fig. The log(/S'e,5o) and 
log(S'e,9o) levels are tabulated in Table ||. We find 
5*6,90 = 5.1 X 10^^ Lq kpc~^ for the radio sample, 
the highest of the three samples. 

In order to determine the effects of the sample 
editing we calculated Se,5o and 5e,9o at the inter- 
mediate editing steps. The most drastic edit- 
ing step was the final one, paring the sample 



down from 244 to 38 sources. The sample be- 
fore this step has log(5e,9o) = 11-29, very close 
to that of the UV sample, but this is due to a 
large number of sources with Se ^ lO^L0kpc-2. 
These are largely removed with the Hll spec- 
tra selection, but also can be removed by apply- 
ing an infrared color selection as done for sam- 
ple II; if no HlI selection is done, and instead 
only galaxies with /eo/Zioo > 0.5 are selected, 
the resulting sample is then very similar to our 
adopted radio sample having log(5e,5o) = 10.09 
and log(S'e,9o) = 11.64. The color selection effi- 
ciently removes cool low surface brightness emis- 
sion from the radio sample, but does not discrim- 
inate against AGN, which like intense starbursts 
tend to have high color temperatures. Although 
we prefer to exclude AGN as best we can, ap- 
plication of the q editing only lowers Se^go by 
0.04 dex. This indicates that allowing radio loud 
AGN to contaminate the sample does not appear 
to significantly effect the sample Se,5o and S'e^go- 

6. Selection effects and biases 

All samples are affected by a surface bright- 
ness selection. The deepest of the UV subsam- 
ples is I.e (Lowenthal et al., 1997). Its limiting 
magnitude for spectroscopy spread over the typi- 
cal Keck FWHM = 0.6" seeing disk, corresponds 
to a limiting ~ 1.4 x lO^LQkpc"^ at z = 3, 
assuming the median value of f3. Samples II and 
III are drawn from the BGS which has a flux 
limit of /so = 5.4 Jy, and is composed primar- 
ily of sources unresolved by IRAS thus having 
FWHM < 2'. For a typical /eo/Zioo = 0.55, this 
corresponds to S'e ~ 3 x lO''L0kpc~^. The re- 
sult that none of the galaxies in samples II and 
III have Se this low is due to the infrared color 
selection (sample II) and Hll spectrum selection 
(sample III). The important point is that sample 
I can not reach the same Se depth as observed in 
the other two samples. 

The biggest systematic flux uncertainty for the 
UV sample comes from how we choose to model 
fig. ||. As noted in §3.1, we expect a ~ 0.4 dex 
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in -Ffir/-^220 from the scatter in Fig. |T[ Indeed 
if we had used the Kinney et al. (1994) extinc- 
tion curve in our analysis we would have derived 
'S'e.sO) 5'e,9o to be lower by 0.2 and 0.4 dex respec- 
tively for sample I. We have not corrected sam- 
ples II and III for the fraction of starburst light 
that escapes dust reprocessing. This should be 
a small correction since the median Li^/Lb ~ 5 
for the galaxies in sample II (Lehnert &: Heck- 
man 1995; Armus et al. 1987). Since much of 
the B band flux in these galaxies probably arises 
outside of the starburst we can expect Se of the 
starburst to be systematically underestimated by 
^ 20%. 

Quantifying the Re biases is more difficult, but 
we can get an idea of the direction they will take. 
Relative to our stated aim of measuring the Re of 
the ionizing population, most of our techniques 
overestimate Re due to the effects of the ISM. In 
principle i?e(UV) provides the most direct esti- 
mate of the stellar population size. It may be in- 
flated by extended nebular continuum and scat- 
tered light. However M95 show that the former 
is probably small (~ 10% in flux). The resolu- 
tion of the nearest starbursts into stars argues 
against scattered light dominating in the UV. 
Differential extinction (e.g. dust mixed in with 
the stars) also may inflate iie(UV) and Re{Roi). 



Although in §3.1 we argue that the UV extinc- 
tion is dominated by a foreground screen (hence 
extinction will not be differential), M95 note that 
an idealized uniform screen can not fully describe 
the dust distribution of most starbursts. ReiJ^a) 
(sample II) may be more strongly effected by dif- 
ferential extinction than i?e(UV) (sample I), de- 
spite the sense of the wavelength difference, be- 
cause the IRAS galaxies are selected to be dusty. 
In addition, sample II. g is selected by edge-on ap- 
pearance. Hence, Ha emission may be extincted 
by cool outer disk dust, far from the starburst, 
resulting in a large i?e(Ha) bias. We may also ex- 
pect that i?e(HQ) to be inflated in galaxies dis- 
playing galactic winds. However it is not clear 
that this is a large effect. On the one hand. 



Marlowe (1997) find that the effect can be a 
large (factor of five in Re) in windy blue compact 
dwarves. On the other hand, using Lehnert's &: 
Heckman's (1995) statistic log(RHa/Ri?) we find 
no noticeable increase of ^e^go for subsample II. g 
when the ten galaxies with the strongest winds 
are excluded from the subsample. As for sam- 
ple HI, the i?e(21cm) measurements will be least 
effected by extinction. A bias towards high sur- 
face brightness will result from using the Gaus- 
sian HWHM to approximate i?e(21cm), instead 
of doing a full curve of growth analysis, and by 
using radio interferometric measurements pref- 
erentially of high angular resolution. Both pro- 
cedures discard low surface brightness extended 
emission. On the other hand, radio emission 
tends to be more extended than star formation 
in nearby galaxies (e.g. Marsh &; Helou, 1995), 
probably due to the propagation of cosmic rays 
in the hosts. 

In summary, selection effects and measure- 
ment biases may significantly effect the Se distri- 
butions. i?e(UV) and Re{^oi) are both likely to 
be inflated with respect to the Re of the ionizing 
population; the latter is likely to be more severely 
affected. Competing physical and measurement 
biases may cause -Re (21cm) to be either under 
or overestimated. We expect systematic biases 
up to a factor of three in UV flux. Agreement 
of 5e,9o values to within this factor indicate con- 
sistent starburst intensity limits as best we can 
determine. 

7. Discussion 

7.1. The Surface Brightness Limit 

Table |^ shows that the (full) samples II and 
HI have the largest difference in S'e^go, 0.8 dex. 
Since they were originally drawn from the same 
source, the IRAS BGS, the factor of six differ- 
ence in 5*6,90 is illustrative of the effect that de- 
tails of final sample selection and measurement 
bias can have. Since the fluxes used to calculate 
Se are essentially the same (tied to Ffir), the 



12 



differences in 5e,9o probably arise in the Re es- 
timates, per flaps for the reasons discussed in 
The 5e,9o of sample I splits the difference of the 
other two samples, agreeing with each to within 
a factor of three. This is a reasonable level of 
agreement as discussed above. 

We conclude that to better than an order 
of magnitude, the maximum surface brightness 
5'e,9o of starbursts are the same, when deduced 
from radio, optical and UV observations. We find 
a mean ^e^go ~ 2.0 x 10^^ Lq kpc~^ (with a fac- 
tor of three uncertainty) as the characteristic sur- 
face brightness of starbursts, found by averaging 
in the log the results of the three samples. The 
lack of a simple wavelength dependence suggests 
that this is a star formation intensity limit, not 
an opacity effect. Starbursts orders of magnitude 
more intense than deduced from UV observations 
are not common. 

Further evidence of the physical significance 
of 5e,9o is given by the central pressure Pq of 
windy FIRGs. For a sample of 12 galaxies 
(mostly starbursts) undergoing strong galactic 
wind Heckman et al. (1990) estimate Pq us- 
ing the [S II]AA6716/6731 A emission line ra- 
tio. They find the mean (it dispersion) Pq = 
2.8 ± 1.2 X 10"^ dy cm . In the starburst mod- 
els of Chevalier & Clegg (1985), Pq of a free flow- 
ing wind is given by Pq = O.llSpR^^, where P IS 
the momentum flux from the starburst, and R 
is the starburst radius, and the constant is for 
cgs units. For our standard starburst population 
model (§ ^.2D the momentum flux due to super- 
novae and stellar winds is given by 



p = 2.1 X 10 



23 



L 



bol r -li 
erg cm J 



(14) 



The constant varies by a factor of 1.6 for burst 
durations in the range of 5 to 100 Myr. The 
star formation intensity required to produce Pq 
is then 



Se 



lOii Lq kpc-2 1.63 X 10-9 dy cm"2 ' 



(15) 



Hence the mean value of Pq measured by Heck- 
man et al. (1990) corresponds to Se = 1.7 x 
lO^^L0kpc^. This is very close to the adopted 
56,90, but derived using methods completely in- 
dependent of those in §^ - §|5|. 

Since we are measuring the surface brightness 
of star formation tracers, we can convert Se,9Q 
to an equivalent star formation intensity {A4 per 
unit area) of 



^e,90 



45A^0kpc ^yr ^ 



(16) 



using eq. ^ (and bearing in mind the caveats of 
Sec. |3.2| ). Strictly speaking the observed emission 
is dominated by high-mass star formation. If the 
lower mass limit of our adopted Salpeter IMF is 
5 Mq, then Se,9o is reduced by a factor of 5.5. 

Because selection effects preferentially depop- 
ulate low Se sources, 90% or more of starbursts 
have <; 2.0 x 10^^ Lokpc"^. Thus it is close 
to an upper limit to the surface brightness of 
starbursts. This implies that starbursts with 
,90 can only be more luminous by be- 
ing larger. The UV results in particular show 
that the same limit applies over two orders of 
magnitude in Re (four in L) and displays little 
evolution out to z ~ 3. 

The star formation intensity is much higher in 
starbursts relative to that integrated over typi- 
cal galactic disks. This is illustrated in the bot- 
tom panel of Fig. Q which shows the Se distri- 
bution of normal disk galaxies, derived from the 
Ha imaging sample of Ryder & Dopita (1993). 
Here bolometric fluxes due to high-mass star for- 
mation were derived from Ha fluxes using the 
Zanstra method. Re estimates were kindly pro- 
vided by S.D. Ryder, and correspond to the ef- 
fective radius of the elliptical aperture contain- 
ing half of the total Ha emission. Table |l| lists 
/S'e,5o, and 5e,9o for this sample. The typical Se of 
star formation in normal spiral galaxies is about 
three orders of magnitudes fainter than that of 
commonly selected starbursts. Much of the dif- 
ference can be ascribed to the filling factor of star 
formation. The Ha morphology of normal spiral 
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galaxies is characterized by numerous, but widely 
separated Hll regions, whereas starburst galax- 
ies are typically dominated by a single central 
super-bright star forming complex. Individual 
giant Hll regions in normal galaxies may have 
Se approaching that found in starbursts. 

7.2. Starburst Regulation 

What are the regulating mechanisms that lim- 
its the global starburst intensity? The central 
pressure estimates discussed above suggests that 
feedback from high-mass star formation may be 
involved in the regulation (cf. the discussion of 
M82 by Lehnert & Heckman 1996). However, 
it is unclear how this would regulate star for- 
mation intensity. The near equivalence between 
i?e(Ha) and the rotation curve turn-over radius 
Rto found by Lehnert & Heckman (1996) sug- 
gests that the global dynamics of the host in lim- 
iting at least the size of the starburst. Here we 
consider whether the limit is related to the 
stability of inner disks of galaxies. 

Following on the work of Toomre (1964) and 
Quirk (1972), Kennicutt (1989) showed that the 
radial gas surface density profiles (i?) in spiral 
galaxies closely follows, but is somewhat under, 
the surface density at which gas disks would be 
unstable to their own self gravity. Specifically, 
Kennicutt defines the critical density for star for- 
mation as 

where a is the velocity dispersion of the gas, a = 
0.67 is an empirically determined constant, and 



R- 



Hr 



1/2 



+ 4S7^ 



(18) 



is the epicyclic frequency for angular frequency 
n = V{R)/R with V{R) being the rotation 
curve. He found that massive star formation in 
Hll regions seems to be inhibited except where 
Tjg > Sc. The so called "Kennicutt law" has 
proven to be very successful in determining the 



location of H II regions (or lack of them) in a va- 
riety of galaxies including the normal spirals, low 
surface brightness galaxies (van der Hulst et al, 
1993), and HIT galaxies (Taylor et al, 1994). The 
overall success of the Kennicutt law (however cf. 
Ferguson et al., 1996; Meurer et al, 1996) impli- 
cates self gravity in regulating star formation in 
gaseous disks. 

Kennicutt's original calibration of a was done 
for galaxies with star formation in the flat part 
of the rotation curve. What happens in the 
rising portion of the rotation curve where star- 
bursts occur? Assuming a simple linear rise, O 
is constant (solid body rotation), and we have 
Sc = aaQ/{l.Q8G). We assume that Sc, 
since in the Kennicutt scenario, once Tig reaches 
or surpasses Sc, star formation becomes rela- 
tively efficient, thus regulating Tg. This assump- 
tion certainly holds in particular cases, such as 
NGC 3504, where Tg traces Sc well into the solid 
body portion of the rotation curve (Kenney et al, 
1993). On causality grounds, the quickest time- 
scale expected for star formation is the dynamical 
time-scale given by 



td 



' 3tt 
16Gp 



(19) 



(Binney & Tremaine, 1987) where p is the mean 
interior density. Note that, for a cold massless 
disk in the potential of a dominant spherical mass 
distribution, the solid body portion of the rota- 
tion curve corresponds to a constant density core. 
In that case tdyn = ^vril"^. The maximum star 
formation intensity expected for gas at Sc is then 

= 9 X 10"^ aan^ [Mq yr'^ kpc'^] , 



Sc < 



td 



yn 



(20) 

where a is in units of km s ^ , and 17 is in units 
of kms~^ kpc^^. This can be expressed in terms 
of bolometric surface brightness using eq. 0: 



Sc ~ 4.0 X 10^ aan^ [Lq kpc^^]. 



f21) 



using the same units as eq. Figure ^ plots Se 
versus Q for the galaxies in Lehnert & Heckman 
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(1996). There is a weak correlation (correlation 
coefficient R = 0.61) in the sense predicted by 
eq. which is plotted for a = 0.67, a = 15 
kms~^ (the central Hi velocity dispersion of a 
representative normal disk galaxy; Dickey et al, 
1990). More importantly, almost all the galax- 
ies in this sample have Se < Sc as expected. 
We conclude that our results are consistent with 
Kennicutt's threshold model of star formation. 
This implies that starbursts may be regulated 
by the same gravitational instability mechanism 
that applies to the outer disks of galaxies. The 
extraordinary nature of starbursts is not due to 
extraordinary physics. Rather it is due to the 
high densities required to reach the instability 
limit, and the consequent short minimum dy- 
namical time-scales (see also Elmegreen, 1994). 

Although application of the Kennicutt law is 
insightful, it does not explain why the same sur- 
face brightness limit appears to hold over two 
orders of magnitude in size. It suggests that 
the blame should be shifted to a fundamental 
limit on the angular frequency or to the cen- 
tral density since po oc Q'^. But we have no in- 
sight into whether such a limit really exists, or 
why it should exist. It also does not explain the 
efficiency of star formation. We conclude that 
the stability of a gaseous disk against self grav- 
ity probably plays an important role in limiting 
starburst intensity, but in and of itself is not suf- 
ficient to explain the Se limit of starbursts. 

7.3. Galaxy formation 

How does 5*6,90 relate to the global star for- 
mation history of normal galaxies? Can they 
be made in a high intensity starburst? Normal 
galaxies also show preferred surface brightnesses. 
Bright spirals obey the Freeman (1970) law hav- 
ing exponential disks with extrapolated B band 
central surface brightnesses of i3(0)c = 21.65 mag 
arcsec"^, which translates to an effective surface 
brightness of Se = 5.1 x IO'^Lb,© kpc~^. Mc- 
Gaugh et al. (1995) showed that this is actually 
an upper limit to the surface brightness distri- 



bution of disk galaxies when selection effects are 
taken into account. Elliptical (E) galaxies fall 
on the so-called "fundamental plane" in log(-Re)) 
log(5e), and log{a). Using the fundamental plane 
projections of Djorgovski & Davis (1987) we es- 
timate the properties of a representative low- 
luminosity E galaxy to be Re = 1.2 kpc, a = 
180kms'\ and Seirc) = 9.0 x lO^L^^c© kpc-^ 
whereas for a representative high-luminosity E 
galaxy we adopt Re = 12 kpc, a = 350kms~^, 
and Seirc) = 1-4 x lO^L^g,© kpc"^. Here the ob- 
servations are defined in the re passband (Djor- 
govski, 1985). 

An estimate of the minimum time needed to 
build a galaxy (e.g. at the starburst intensity 
limit) is given by 



^build 



{M/L)Se 



Myr, 



(22) 



^e,90 



where Sg^go is given in eq. |T^. Taking M/Lb ~ 
2.5 to be appropriate for spirals (Puche & Carig- 
nan, 1991) and M/L^q ~ 5 for E galaxies (Djor- 
govski & Davis, 1987) yields tbuiid ~ 3, 98, and 
15 Myr for spirals, low luminosity E and high 
luminosity E galaxies respectively. 

As before, the minimum time-scale for star 
formation will be on the order of the dynamical 
time. This is given by 



2Re 



a 



2Re 



(23) 



(following equation 4-80b in Binney & Tremaine, 
1987). Taking Re = 1.68a^^ (Freeman, 1970) 
and using the data of Puche & Carignan (1991) 
typically t^yn ~ 60 Myr for spiral galaxy disks. 
Using the representative E galaxy parameters 
listed above, tdyn ~ 13, 67 Myr for low and high 
luminosity E galaxies respectively. 

We see that tbuiid is on the order of, or larger 
than, tdyn for E galaxies, and so it is plausible 
that they were made by a maximum intensity 
starburst - this scenario would not violate burst 
duration limits set by causality. It is thought 
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that the epoch of E galaxy formation is at 2 ^ 1 
(e.g. Bender et ai, 1996), hence it is plausible 
that the high-z starbursts are elliptical galaxies 
forming. This issue is further addressed below. 

Spiral disks, on the other hand, have tbuiid ^ 
*dyn) so a maximum intensity burst would pro- 
duce many more stars in a dynamical time-scale 
than are observed in spiral disks. These consider- 
ations are in accord with observations of present 
day disk galaxies which indicate that star for- 
mation occurs on a few Gyr or longer time-scale 
(Kennicutt et al, 1994). Although maximum in- 
tensity bursts are too strong to build disks, disk 
galaxies usually have central high surface bright- 
ness bulges. Bulges fall on the fundamental plane 
and thus they can be made by maximum inten- 
sity bursts. Indeed, the surface brightness pro- 
files of central starburst galaxies are somewhat 
akin to those of spiral galaxies, but with the 
"bulge" being blue and corresponding to the star- 
burst (cf. Schade et al. 1996; Lilly et al. 1995). 

7.4. Extinction in high redshift galaxies 

In §^ we showed that there is little apparent z 
evolution in Se of the UV-selected sample. One 
of the most important corrections we make to the 
UV data is the extinction correction which turns 
out to be substantial. Here we explore whether 
these corrections are reasonable, and the impli- 
cations arising from the extinction correction. 

Figure ^ shows the distribution of /? values for 
the low and high redshift galaxies. Here we plot 
only one point in the cases where a galaxy was 
separated into multiple starburst components, 
since the same f3 estimate was used for all compo- 
nents in these cases. The arrows show the intrin- 
sic Po = —2.5 for an unreddened ionizing pop- 
ulation, and the median /? of the sample. The 
median shifts by less than 0.05 in /? if the local 
sample is excluded. Furthermore, the range of 
observed /? is virtually identical at high and low 
redshifts. The high- 2; starbursts thus have very 
similar UV color distributions to the low-z sam- 
ple. 



Figure illustrates that most starbursts, in- 
cluding the high- 2: ones, are redder than naked 
ionizing populations. The high-z samples are se- 
lected to have the bluest possible rest-frame UV 
colors. Young unreddened populations will pref- 
erentially be selected if they exist. Either there 
are very few high-z galaxies dominated by ion- 
izing populations, and/or some of the ionizing 
populations are reddened. 

Opacity from the Lya forest may redden the 
photometric (3 values, if the blue filter in the 
color index samples rest wavelengths bluewards 
of 1216A. This is likely to be a problem only for 
subsample I.e. Excluding this subsample changes 
the median /3 by less than 0.05. Hence Lya forest 
opacity is not the reddening mechanism. 

Our interpretation is that the reddening is due 
to dust extinction. We know this is valid for the 
local UV sample (see Fig. |l|). Unfortunately the 
confirmation spectra of the high-z galaxies do not 
extend to the rest-frame optical where the emis- 
sion lines should be strong, allowing confirma- 
tion of the presence of an ionizing population. 
The same rest frame UV absorption features are 
seen in both the low and high-z samples, how- 
ever they are relatively narrow in the high-z sam- 
ple (Lowenthal et al, 1997), so it is not clear 
whether they have a stellar or interstellar origin. 
Hence, we can not rule out the possibility that 
the high-2 sample is contaminated by older non- 
ionizing stellar populations. It is unlikely that 
there is significant contamination by faded star- 
bursts, because aging has a more severe effect 
on luminosity than our extinction corrections. It 
will take ~ 400 Myr after a starburst turns off for 
it to redden to the observed median /3 = —1.1, 
during that time it will have faded by a factor 
of ~ 500 (assuming the burst duration <^ 400 
Myr; see Fig. 16 of M95), whereas the extinc- 
tion correction for this /? is a factor of eight in 
the F220W band. Hence, post-bursts would have 
been that much brighter when they were ionizing 
populations, and such ultra-luminous starbursts 
are not detected at high-z. More worrisome is 
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contamination by extended duration star form- 
ing populations. However, durations longer than 
300 Myr are required to reach the median (3 even 
if the IMF slope is a = 3.3 (LH95). If contami- 
nation is severe we would have the curious effect 
that S'e and (3 evolution "conspire" to mimic no 
evolution and a simple extinction correction. We 
prefer the simpler interpretation that the same 
physical mechanism, extinction by dust, reddens 
both the low and high-z starbursts. 

The top axis of Fig. |^ translates (5 to extinc- 
tion at 1620A (the approximate rest wavelength 
of the HDF observations longward of the Lyman 
Limit) using the C94 extinction law. The median 
extinction is 2.9 magnitudes at this wavelength (a 
factor of 15 in flux). For the Kinney et al. (1994) 
extinction law, the median extinction would be 

2.0 magnitudes. In either case, the extinction is 
substantial. 

In their pioneering paper, Madau et al. (1996) 
use these high-z galaxies to evaluate the metal 
enrichment history of the early universe deriving 
apparent metal ejection rates pz{z) = 3.6, 0.62, 

1.1 xlO-'^Moyr-iMpc-i at z = 2.75,3.25,4 
(see Madau et al. for details). They note that 
these values are underestimates because (1) they 
do not sample the full luminosity function, (2) 
they represent only the least dusty systems, and 
(3) they assume no dust extinction in their anal- 
ysis. We can improve on their work by cor- 
recting their results for the median ^i620 = 2.9 
mag extinction seen in starbursts. This yields 
pz{z) = 54, 9.4, 17 X 10"^ M© yr'^ Mpc'^ for the 
above three redshifts. Strictly speaking these re- 
sults are still lower limits to the true pz{z) of the 
universe. 

Figure ^ shows our revised version of the pz{z) 
evolution diagram (cf. Fig. 9 of Madau et al.). 
Here we adopt the flat cosmology used by Madau 
et al. Hq = 50 kms~^ Mpc~^, go = 0.5, and as- 
sume A = 0. The X axis here is look-back time in 
linear units, and the y axis is linear in pz{z). All 
data points are shown as lower limits: the high-z 
data because completeness corrections still have 



not been addressed; and the data with z < 1 be- 
cause they have not been corrected for extinction. 
If the extinction correction to the other points is 
small, then our new estimates suggest that the 
epoch of peak galaxy formation may be as early 
as z ~ 3. Clearly much work needs to be done to 
put all estimates of pz{z) on the same scale. 

Our new estimates of pz{z) are well above the 
mean metal ejection rate of the universe {pz) ~ 
4.2 X lO-^Moyr-^Mpc-i (Madau et al, 1996). 
Strictly speaking it too is a lower limit, since it 
neglects low surface brightness galaxies and in- 
tergalactic and intracluster gas in the estimation. 
Mushotzky & Loewenstein (1997) provide an es- 
timate of {pz) = 3.4 X lO-^Mgyr-^Mpc"^ in 
the redshift interval z = 1 to 6 - the epoch of 
elliptical galaxy formation. This level is plotted 
as a dashed line in Fig. |^. In Madau et al. it 
seemed possible that the universe was relatively 
quiescent at high redshifts. It now appears that 
the universe was forming stars, and producing 
metals, at a high rate at z > 2. Our results 
imply that rather than a quiescent universe we 
are observing an active but moderately obscured 
early universe. 

8. Summary 

We have examined a diverse set of observa- 
tions of starburst galaxies, spanning six decades 
of wavelength, sizes ranging from ~0.1 to 10 kpc, 
and redshifts out to z = 3.5. The observations 
were used to estimate effective radii Re, total lu- 
minosities L, and effective surface brightnesses 
Se on the bolometric scale. The samples were 
drawn from observations obtained in 1. the rest 
frame space ultraviolet (UV), 2. the far- infrared 
and Ha (FIR/Ha), and 3. 21cm radio emission. 
The UV data, which spans the largest range in 
z, were consistently corrected for dust extinction 
as estimated from the UV reddening of the star- 
bursts. We have demonstrated that this correc- 
tion adequately models the redistribution of ra- 
diation from the UV to the FIR in nearby star- 
bursts. 
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The extinction in the UV is significant. High- 
z galaxies observed in tlic rest frame UV are 
too red to be naked starbursts. If this red- 
dening is due to dust, the typical implied ex- 
tinction is ~ 2 to 3 mag at A = 1620A. The 
corrected metal production rate at z ^ 3 is 
then ^ 2 X 10~^ Mpc~^ yr~^, about five 
times or more higher than Madau et al.'s (1996) 
estimate of the HTibble time averaged {pz) ~ 
4.2 X 10^"^ A^0 Mpc^'^yr^^, but in good agree- 
ment with Mushotzky &: Loewenstein's (1997) 
estimate of (pz) ~ 3.4 x 10~^ Mq Mpc"^ yr~^ 
during the epoch of elliptical galaxy formation. 
Depending on the amount of extinction required 
for the data having z < 2, the peak in pz{z) 
could be pushed back to z ^ 2. Applying plausi- 
ble extinction corrections illustrates that rather 
than being quiescent, much of the star formation 
in the early universe may be obscured. 

The upper limit to the surface brightness of 
starbursts is parameterized by Se,9o - the 90th 
percentile surface brightness of the samples. This 
statistic is susceptible to systematic uncertain- 
ties due to measurement biases and the heteroge- 
neous sample selection. We find that the samples 
have the same 5*6,90 to within a factor of three. 
This is about the level of agreement we expect for 
intrinsically similar samples affected by known 
flux biases. We conclude that to within a fac- 
tor of a few the same S'e.go applies to both dusty 
and relatively unobscured starbursts, it does not 
vary with size for i?e ^ 0.1 kpc, and shows little 
evolution out to z 3.5. 

Our adopted 90th percentile surface bright- 
ness limit is the mean value (in the log) ^g^go ~ 
2.0 X 10^^ Lq kpc~^ of the three samples. The 
absence of a strong wavelength variation of Se,9o 
indicates that it corresponds to a physical limit 
on star formation intensity, rather than being an 
opacity effect. Further evidence of the physi- 
cal significance of 5e,9o comes from the central 
ISM pressure Pq = 2.8 ± 1.2 x lO^^^dycm-^ 
measured in galaxies (mostly starburst) under- 
going a strong galactic wind by Heckman et al. 



(1990). Using the models of Chevalier &: Clegg 
(1985) and Leitherer Sz Heckman (1995), the 
corresponding surface brightness of the stellar 
population pressurizing the ISM is Sg = 1.7 x 
lO^^L0kpc^. This is identical to S'e.go, but is 
derived from completely different considerations. 

Our adopted Se,9o corresponds to a star for- 
mation intensity of ile,go ~ 45A^0 kpc~^ yr~^. 
Since starbursts represent the most intense star 
formation in the universe, and since selection 
effects depopulate low Se sources, these values 
of S'e.go and Se,go parameterize the maximum 
global intensity of star formation observed in the 
universe. Within starbursts there are localized 
sources with Se > 2 X lOiiL0kpc-2. These 
are star clusters. However, their integrated light 
does not dominate starburst emission, and their 
sizes are limited to Re <^ 10 pc. The intense clus- 
ter forming mode of star formation apparently 
does not operate for i?e ^ 0.1 kpc. 

Our results imply that some robust mecha- 
nism is limiting the global intensity of starbursts. 
The Pq comparison implicates mechanical energy 
feedback from supernovae and stellar ejecta in 
regulating starbursts. We also examined the con- 
sequences of disk stability models (e.g. Kenni- 
cutt, 1989) to the central portions of galaxies 
where the rotation curve is solid body like and 
most starbursts are observed to occur. There 
we expect an upper limit to Se oc aU^ oc ap, 
where a is the gas velocity dispersion in the disk, 
ri is the angular frequency, and p is the spher- 
ical volume density. We show that most star- 
bursts with rotation curve data obey this limit 
for an adopted central velocity dispersion o" = 15 
kms~^. This model implies that starbursts have 
higher Se than normal disk galaxies because the 
star formation occurs in the center of galaxies 
where the critical density is high and the star 
formation time-scale is short. Unfortunately, this 
model (while illuminating) does not explain the 
constancy of S'e.go, it only casts it in terms of a 
limit on Q or the central value of p. 

Normal galactic disks could not have been 
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made in a maximum intensity burst because it 
would have been accomplished in much shorter 
than the dynamical time scale tdym the mini- 
mum time set by causality arguments. How- 
ever, bulges, which often dominate the center of 
disk galaxies, can be made in a maximum inten- 
sity burst without violating such causality limits. 
The same goes for elliptical galaxies. Consid- 
ering the Re and L of the high-z galaxies, and 
the agreement of the corrected metal production 
rate derived from them with that expected at 
the epoch of E galaxy formation, it is tempting 
to speculate that at z ~ 3 we are witnessing the 
formation of elliptical galaxies (Giavalisco et al. 
1995). 

Acknowledgements: Numerous astronomers watched 
this work develop and provided useful hints of 
samples to include, and other suggestions. In 
particular we would like to thank Mark Dick- 
inson, Mauro Giavalisco, Jeff Goldader, Piero 
Madau, Stuart Ryder, Nick Scoville, Min Yun, 
and Esther Zirbel. Asao Habe, Chris Mihos, 
Colin Norman, and Rosie Wyse provided useful 
comments on earlier drafts of this paper. We are 
grateful to Alexei Filippenko for obtaining the 
QNY1:32 spectrum for us. Jim Condon saved 
us a lot of effort by emailing us data from C90. 
We are grateful for the support we received from 
NASA through HST grant number GO-05491.01- 
93A and LTSA grant number NAGW-3138. This 
paper is in large part based on observations with 
the NASA/ESA Hubble Space Telescope obtained 
at the Space Telescope Science Institute, which 
is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. Literature searches were 
performed using NED, the NASA/IPAC Extra- 
galactic Database, a facility operated by the Jet 
Propulsion Laboratory, Caltech, under contract 
with the National Aeronautics and Space Admin- 
istration. 



REFERENCES 

Allen, C.W. 1973, Astrophysical Quantities, Third 
Edition, (Athlone, London) 

Armus, L., Hcckman, T.M., & Miley, G.K. 1987, 

AJ, 94, 831 

Armus, L., Heckman, T.M., & Miley, G.K. 1990, 

ApJ, 364, 471 

Bender, R., Ziegler, B., & Bruzual, G. 1996, ApJ, 
463, L51 

Binney, J., & Tremaine, S. 1987, Galactic Dy- 
namics, (Princeton University Press, Prince- 
ton) 

Bower, G.A., Wilson, A.S 1995, ApJS, 99, 543 

Boyle, B.J., Fong, R., Shanks, T., and Peterson, 
B.A 1990, MNRAS, 243, 1 

Burstein, D., & Hciles, C. 1984, ApJS, 54, 33 

Calzetti, D., Kinney, A.L., & Storchi-Bergmann, 
T. 1994, ApJ, 429, 582 (C94) 

Calzetti, D., Kinney, A.L., & Storchi-Bergmann, 
T. 1996, ApJ, 458, 132 

Carozzi, N. 1977, A&A, 55, 261 

Carozzi-Meyssonnier, N. 1978, A&A, 63, 415 

Chevalier, R.A., & Clegg, A.W. 1985, Nature, 
317, 44 

Condon, J.J., Helou, G., Sanders, D.B., & Soifer, 
B.T. 1990, ApJS, 73, 359 (C90) 

Condon, J.J., Huaung, Z.-P. Yin, Q.F., &; Thuan, 
T.X. 1991, ApJ, 378, 65. 

Conti, P.S., Vacca, W.D. 1994, ApJ, 423, L97 

Dickey, J.M., Hanson, M.M., & Helou, G. 1990, 
ApJ, 352, 522 

Disney, M.J. 1976, Nature, 263, 573 

Djorgovski, S. 1985, PASP, 97, 1119 

Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 

Durret, F. & Bergeron, J. 1988, A&AS, 75, 273 

Elmegreen, B.C. 1994, ApJ, 425, L73 



19 



Ferguson, A.M.N., Wysc, R.F.G., Gallagher, J.S. 
Ill, Hunter, D.A. 1996, AJ, 111, 2265 

Freeman, K.C. 1970, ApJ, 160, 811 

Gallego, J., Zamorano, J., Aragon-Salamanca, 
A., & Rego, M. 1995, ApJ, 455, LI 

Giavalisco, M., Maccheto, F.D., Madau, P., & 
Sparks, W.B. 1995, ApJ, 441, L13 

Giavalisco, M., Steidel, C.C., & Macchetto, F.D. 
1996, ApJ, 470, 189 

Heckman, T.M., Armus, L., & Miley, G.K. 1990, 
ApJS, 74, 833 

Helou, G., Soifer, B.T., & Rowan-Robinson, M. 

1985, ApJ, 298, L7 

Helou, G., Khan, I.R., Malek, L., &: Boehmer, L. 
1988, ApJS, 68, 151 

Helou, G. 1991, in The Interpretation Of Mod- 
ern Synthesis Observations Of Spiral Galax- 
ies, edited by N.C. Duric, & P. Crane, (Astro- 
nomical Society of the Pacific, San Francisco), 
p. 125 

Holtzman, J. A., Faber, S.M., Shaya, E.J., Lauer, 
T.R., Groth, E.J, Hunter, D.A., Baum, W.A., 
Ewald, S.P., Hester, J.J., Light, R.M., Lynds, 
R., O'Neil, E.J. Jr., Westphal, J.A. 1992, AJ, 
103, 691 

Holtzman, J.A., Watson, A.M., Mould, J.R., 
Gallagher, J.S. Ill, Ballester, G.E., Burrows, 
C.J., Clarke, J.T., Crisp, D., Evans, R.W., 
Griffiths, R.E., Hester, J.J., Hoessel, J.G., 
Scowen, P.A., Stapelfeldt, K.R., Trauger, J.T., 
Westphal, J.A. 1996, AJ, 112, 416 

Hunter, D.A., O'Connell, R.W., Gallagher, J.S. 

1994, AJ, 108, 84 

IRAS Point Source Catalog, ver. 2. 1988. Wash- 
ington: GPO 

Kenney, J.D.P., Carlstrom, J.E., &; Young, J.S. 

1993, ApJ, 418, 687 

Kennicutt, R.C. 1989, ApJ, 344, 685 

Kennicutt, R.C., Tamblyn, P., & Congdon, C.E. 

1994, ApJ, 435, 22 



Kim, D.-C, Sanders, D.B., Veilleux, S., Maz- 
zarella, J.M., k Soifer, B.T. 1995, ApJS, 98, 
129 

Kinney, A.L., Calzetti, D., Bica, E., &; Storchi- 
Bergmann, T. 1994, ApJ, 429, 172 

Lehnert, M., & Heckman, T.M. 1995, ApJS, 97, 
89 

Lehnert, M., & Heckman, T.M. 1996, ApJ, 472, 

546 

Leitherer, C, & Heckman, T.M. 1995, ApJS, 96, 
9 (LH95) 

Leitherer, C, Ferguson, H.C., Heckman, T.M., 
Lowenthal, J.D. 1996, ApJ, 454, L19 

Lilly, S.J., Hammer, F., Le Fevre, O., & Cramp- 
ton, D. 1995, ApJ, 455, 75 

Lilly, S.J., Le Fevre, O., Hammer, F., &; Cramp- 
ton, D. 1996, ApJ, 460, LI 

Lisenfeld, U., Volk, H.J., k Xu, C. 1996, A&A, 
306, 677 

Lowenthal, J.D., Koo, D.C., Guzman, R., Gal- 
lego, J., Phillips, A.C., Faber, S.M., Vogt, 
N.P., Illingworth, G.D., & GronwaU, C. 1997, 
ApJ, in press 

Madau, P., Ferguson, H.C., Dickinson, M.E., Gi- 
avalisco, M., Steidel, C.C., & Fruchter, A. 
1996, MNRAS, 283, 1388 

Maoz, D., Barth, A.J., Sternberg, A., Filippcnko, 
A.V., Ho, L.C., Macchetto, F.D., Rix, H.-W., 
Schneider, D.P. 1996a, AJ, 111, 2248 

Maoz, D., Filippenko, A.V., Ho, L.C., Mac- 
chetto, F.D., Rix, H.-W., Schneider, D.P. 
1996b, ApJS, 107, 215 

Marlowe, A.T. 1997, Ph.D. Thesis, The Johns 

Hopkins University 

Marsh, K.A., & Helou, G. 1995, ApJ, 445, 599 

Mathis, J.S. 1990, ARA&A, 28, 37 

McGaugh, S.S., Bothun, G.D., & Schombert, 
J.M. 1995, AJ, 110, 573 

Meurer, G.R., Freeman, K.C, Dopita, M.A., & 
Cacciari, C. 1992, AJ, 103, 60 



20 



Meurcr, G.R., Heckman, T.M., Leitherer, C, 
Kinney, A., Robert, C, & Garnett D.R. 1995, 
Astron. J., 110, 2665 (M95) 

Meurer, G.R., Carignan, C, Beaulieu, S., & 
Freeman, K.C. 1996, AJ, 111, 1551 

Mushotzky, R.F., & Loewenstein, M. 1997, ApJ, 
accepted 

O'Connell, R.W., Gallagher, J.S. Ill, Hunter, 
D.A., Colley, W.N. 1995, ApJ, 446, LI 

O'Connell, R.W., & Mangano, J.J. 1978, ApJ, 
221, 62 

Peterson, C.J. 1980, PASP, 92, 397 

Puche, D., & Carignan, C. 1991, ApJ, 378, 487 

Quirk, W.J. 1972, ApJ, 176, L9 

Rubin, V.C., Thonnard, N., Ford, W.K., & 
Burstein, D. 1982, ApJ, 261, 439 

Ryder, S.D., & Dopita, M.A. 1993, ApJS, 88, 415 

Salpeter, E.E. 1955, ApJ, 121, 161 

Sanders, D.B. & Mirabel, I.F. 1996, ARA&A, 34, 
749 

Schade, D., Lilly, S.J., Lc Fcvrc, O., Hammer, 
F., & Crampton, D. 1996, ApJ, 464, 79 

Schweizer, F., Miller, B.W., Whitmore, B.C., 
Fall, S.M. 1996, AJ, 112, 1839 

Seaton, M.J. 1979, MNRAS, 187, 73p 

Soifer, B.T., Sanders, D.B., Madore, B.F., Neuge- 

bauer, G., Daniclson, G.E., Elias, J.H., Lons- 
dale, C.J., & Rice, W.L. 1987, ApJ, 320, 238 

Soifer, B.T, Boehmer, L., Neugebauer, G., & 
Sanders, D.B. 1989, AJ, 98, 766 

Stark, A.A., Gammie, C.F., Wilson, R.W., Bally, 
J., Linke, R.A., Heiles, C, Hurwitz, M. 1992, 
ApJS, 79, 77 

Stcidcl, C.C. & Hamilton, D. 1992, AJ, 104, 941 

Stcidcl, C.C, Pettini, M., & Hamilton, D. 1995, 
AJ, 110, 2519 

Steidel, C.C, Giavalisco, M., Pettini, M., Dick- 
inson, M., Adelberger, K.L. 1996a, ApJ, 462, 
L17 



Steidel, C.C, Giavalisco, M., Dickinson, M., 
Adelberger, K.L. 1996b, AJ, 112, 352 

Taylor, C.L., Brinks, E., Pogge, R.W., & Skill- 
man, E.D. 1994, AJ, 107, 971 

Toomrc, A. 1964, ApJ, 139, 1217 

van der Hulst, J.M., Skillman, E.D., Bothun, 
G.D., McGaugh, S.S., & de Blok, W.J.G. 
1993, AJ, 106, 548 

Veilleux, S., Kim, D.-C, Sanders, D.B., Maz- 
zarclla, J.M., & Soifer, B.T. 1995, ApJS, 98, 
171 

Watson, A.M., Gallagher, J.S., III, Holtzman, 
J.A., Hester, J.J., Mould, J.R., Ballaster, G.E. 
1996, AJ, 112, 534 

Whitmore, B.C., Schweizer, F., Leitherer, C, 
Borne, K., Robert, C. 1993, AJ, 106, 1354 

Whitmore, B.C., & Schweizer, F. 1995, AJ, 109, 

960 

Williams, R.E., Blacker, B., Dickinson, M., Dixon, 
W.V.D., Ferguson, H.C, Fruchter, A.S., Gi- 
avalisco, M., Gilhland, R.L., Hcycr, I., Kat- 
sanis, R., Levay, Z., Lucas, R.A., McElroy, 
D.B., Petro, L., Postman, M., Adorf, H.-M., 
& Hook, R.N. 1996, AJ, 112, 1335 

Zhang, X., Wright, M., & Alexander, P. 1993, 
ApJ, 418, 100 



This 2-column preprint was prepared with the AAS IM^jX 
macros v4.0. 



21 



Table 1 

foc observations of moderate redsiiift starbursts 



Galaxy 


z 


W-F342W 




A- 




alh 


D 


log(Lboi/i0) 


Re 






(STMAG) 


(mag) 


(mag) 


(") 




(Mpc) 




(kpc) 


QNY1:32 


0.4424 


19.93 


0.16 


2.2 


0.32 


2.03 


1338 


12.06 


1.70 


SGP1:10 


0.4237 


19.87 


0.17 


0.0 


0.20 


1.70 


1294 


11.18 


1.14 



^Galactic extinction evaluated at observed A = 3403A. 
^'Intrinsic extinction evaluated at rest Aq = 2320A. 



Table 1 

Statistics of samples and subsamples 



(sub)sample 


N 


log(S^i„)^ 


log(5e,50)" 


log(Se,90)"'^ 


log(.Smax)^ 


Description 


I. 


42 


9.31 


10.49 


11.31 


11.94 


Total UV sample 


La 


11 


9.47 


10.53 


11.32 


11.94 


Local UV sample (M95) 


Lb 


2 


9.93 


10.10 




10.26 


New moderate z starbursts 


I.c 


6 


9.75 


10.23 




11.39 


High-2, Steidel et al. (1996a) 


I.d 


8 


10.11 


10.70 




11.46 


HDF, Steidel et al. (1996b) 


I.e 


15 


9.31 


10.12 


10.63 


10.89 


HDF, Lowenthal et al. (1997) 


I.c-e 


29 


9.31 


10.43 


11.21 


11.40 


All high- 2; UV subsamples 


II. 


48 


8.23 


10.24 


10.92 


11.50 


FIR/Hq sample 


Il.f 


18 


9.83 


10.52 


10.08 


11.51 


Armus et al. (1990) 


Il.g 


30 


8.23 


9.69 


10.57 


10.90 


Lehnert & Heckman (1995,1996) 


III. 


38 


8.55 


10.33 


11.71 


12.32 


21cm obs. Condon et al. (1990) 


IV. 


79 


10.27 


> 11.92'= 


> 12.73^ 


> 13.82'= 


Star clusters in the UV (M95) 


V. 


34 


7.37 


7.22 


7.69 


8.22 


Normal disks: Ryder & Dopita (1993) 



^Se units are Lq kpc~ 

'^Calculated only if (sub)sample size is > 10. 
'^lower limits due to large number of Rg upper limits. 
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Fig. 1. — The ratio of far- infrared to UV fluxes compared to ultraviolet spectral slope /?, adapted from 
Fig. 6 of M95. The data correspond to UV-selected starbursts observed by IRAS and HST or lUE. Only- 
galaxies with isophotal diameters D25 ^ 5' are plotted so as to exclude galaxies with large F220 aperture 
corrections. The hatched region shows the expected /3 for naked ionizing populations. The dashed line 
shows the expected relationship for a starburst having an intrinsic /3o = —2.5 that is reddened and 
extincted by a foreground screen of dust with the Calzetti et al. (1994) extinction law. 
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Fig. 2. — The temporal evolution of F220W luminosity as a fraction of the bolometric luminosity for a 
naked starburst population. The curves were derived from the models of Leitherer and Heckman (1995), 
assuming a Salpeter (1955) IMF slope, between mass limits of 1 and 100 Mq. The effect of extending the 
mass range down to 0.1 Mq is negligible. The solid line is for a constant star formation rate population, 
of duration t, and the dotted line is for an instantaneous starburst with age t. 
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Fig. 3. — Bolometric luminosity Lboi and effective surface brightness Sg plotted against effective radius 
Re for UV-selected starbursts. The dotted and dashed lines correspond to Sg^do and Se.50 of the combined 
sample. The correspondence between symbols and subsamplcs is as follows: circles - La. local starbursts 
(M95); triangles - Lb. moderate-z starbursts (this work); six pointed stars - I.e. high-z starbursts (Gi- 
avalisco et al. 1996; Steidel et al, 1996a); five pointed stars - Ld. HDF high-z galaxies Steidel et al. 
(1996b) sample; four pointed stars - Le. HDF high-z galaxies DEEP sample (Lowenthal et al, 1997). 
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Fig. 4. — Bolometric surface brightness distributions of the UV, FIR/Ha, and radio samples compared 
to that derived from Ha observations of normal galaxies (Ryder &; Dopita, 1993). 
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Fig. 5. — Same as Fig. 3 for star clusters within starbursts (from M95). Arrows indicate upper limits to 
Re- The dotted and dashed lines have the same position as in Fig. 3. 
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Fig. 6. — Same as as Fig. 3 for the FIR/Ha sample, with the dotted and dashed hnes showing 56,90 
and Se,5o hmits of this sample. Symbol - subsample correspondence: squares - II. f. (Armus et al. 1990); 
diamonds - Il.g. (Lehnert &; Heckman, 1995, 1996). 
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Fig. 7. — Bolometric surface brightness versus angular frequency O in the rising portion of the rotation 
curve for galaxies in subsample Il.g (Lehnert k. Heckman 1996). Here inclinations were estimated from 
the axial ratios using cosi h/a. The dotted line shows the expected upper limit expected for gas at 
Tic forming stars over a dynamical time-scale. Since full data are missing for much of subsample Il.g we 
supplemented their kinematics data with published rotation curve data, and have taken R^. = i?to (the 
rotation curve turnover radius) for galaxies with missing Rf, (required to estimate ^e). The supplemental 
data came from the following sources: Armus et al. (1990), Carozzi (1977), Carozzi-Meyssonnier (1978), 
Durret & Bergeron (1988), O'Conneh & Mangano (1978), Peterson (1980), Rubin et a/. (1982), and Zhang 
et al. (1993). 
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Fig. 8. — Distribution of UV spectral slopes (3 for the UV sample. The symbols are the same as for 
Fig. 3. Only one point is plotted where a galaxy resolves into multiple starburst knots. 
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Fig. 9. — Evolution of metal production rate as a function of look-back time (bottom axis), and redshift 
(top axis). The symbol correspondence is as follows: squares - Madau et al. (1996); x - Steidel et al. 
(1996a); circles - Lilly et al. (1996); triangle - Gallego et al. (1995). The data points are all lower limits. 
For z > 2 this is because completeness corrections have not been made. For z <1 this is because no dust 
extinction corrections have been made. The dotted line shows the Hubble time averaged metal production 
rate estimated by Madau et al. (1996). The dashed line segment shows the mean metal production rate 
during the epoch of elliptical galaxy formation as estimated by Mushotzky & Loewenstein (1997). 
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